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' ^ ' The nucleon form factors of the energy momentum tensor (EMT) [1] were subject to modest interest in literature 
I I for a long time - probably because the only known process, where they (in principle) could directly be "measured" , is 
elastic scattering of gravitons off the nucleon. The situation changed, however, with the advent of generalized parton 
distribution functions (GPDs) [2-5] accessible in hard exclusive reactions [6-13], see [14-18] for reviews. The form 
y—i • factors of the quark part of the EMT of QCD — we use the notation Af^(i), J^{t) and df{t), see the definition below 
^ ' in Eq. (1) — appear as certain Mellin moments of the unpolarizcd quark GPDs. 

' The form factor M^(0) is known at zero- momentum transfer from inclusive deeply inelastic scattering experiments, 
^ID I telling us that quarks carry only about half of the momentum of (a very fast moving) nucleon, and that the rest is 

■ carried by gluons. The appealing perspective is to access by means of GPDs information on J'3(t), which — after 
I extrapolating to zero momentum transfer t = — would reveal how much of the nucleon spin is due to quarks [3] . The 

f — . third form factor, di{t), is equally interesting — promising to provide information on the distribution of strong forces 

■ in the nucleon [19, 20] similarly as the electromagnetic form factors contain information about the electric charge 
distribution [21]. The information content encoded in GPDs is, however, far reacher than that, see Refs. [22-24]. 

' rS i' In this work we study the form factors of the EMT in the framework of the chiral quark soliton model (CQSM) 
I I [25, 26]. The model provides a field theoretic description of the nucleon in the limit of a large number of colours Nc, 
where the nucleon appears as a chiral soliton of a static background pion field [27]. Numerous nucleonic properties, 
(— I I among others form factors [28-31] usual quark and antiquark distribution functions [32-37] and GPDs [38-43], have 
been described in this model without adjustable parameters. As far as those quantities are known an agreement with 
phenomenology was observed typically to within an accuracy of (10 — 30)%. 

Our study provides several new results. In particular, we compute the spatial density distributions, and mean 
^ ' square radii of the operators of different components of the energy momentum tensor and its trace. This provides 
insights not only on, for example, how the "mass" or the "angular momentum" are distributed in the nucleon. Of 
particular interest are the results for the spatial distribution of strong forces in the nucleon. As a byproduct we 
learn how the soliton acquires stability in the GQSM. We also observe a physically appealing connection between the 
criterion for the stability of the nucleon, and the sign of the form factor di{t) at zero-momentum transfer. 

We present results for the form factors which are of practical interest especially in the case of J'^{t). As exclusive 
reactions yield information on J^{t) only at finite < < 0, some guidance from reliable model calculation might be of 
interest for the extrapolation i — > required to conclude how much quarks contribute to the nucleon spin. 

We explore the chiral character of the model to study chiral properties of the form factors. In particular, we derive 
the leading non- analytic chiral contributions to the form factors in the large Nc limit. These non-analytic (in the 
current quark mass) terms are model-independent. In fact, our results coincide with results from chiral perturbation 
theory [44-46], provided one takes into account that the latter is formulated for finite iVc = 3 [47, 48]. 

The implicit pion mass dependence of the form factors is of interest in the context of the chiral extrapolation of 
lattice QCD data [49-54]. This topic can be addressed in the CQSM [55] which will be done in a separate work [56]. 

The note is organized as follows. Sec. II provides a general discussion of the EMT form factors. Sec. Ill introduces 
the model. In Sec. IV we derive the model expressions for the form factors, and discuss the numerical results for 
the densities of the static EMT in Sees. V, VI, VII. In Sec. VIII we present the results for the form factors, and 
conclude our findings in Sec. IX. The Appendices contain: a digression on alternative notations, a discussion of 
general properties of the densities of the static EMT, technical details on the model expressions, and explicit proofs 
for the consistency of the model. 
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II. FORM FACTORS OF THE ENERGY-MOMENTUM TENSOR 



The nucleoli matrix element of the symmetric energy-momentum tensor of QCD is characterized by three scalar form 
factors [1, 3]. The nucleon matrix elements of the quark and gluon parts of the symmetric QCD energy- momentum 
tensor (EMT) can be parameterized as [3, 19] (see App. A for an alternative notation) 



{p'\T^fm\p) = Hp') 



u 



ip)- (1) 



Here {T^u) t'^'^ quark (gluon) part of the QCD energy-momentum tensor. The nucleon states and spinors are 
normalized by {p'\p) = 2p^ {2ttY 5^^\y)' — p) and u(p)u{p) = 2Mjq where we suppress spin indices for brevity. The 
kinematical variables are defined as P = {p + p')/2, A = {p' — p), t = A^. The form factor c{t) accounts for non- 
conservation of the separate quark and gluon parts of the EMT, and enters the quark and gluon parts with opposite 
signs such that the total (quark+gluon) EMT is conserved. 

The nucleon form factors of the EMT arc related to the unpolarizcd GPDs iJ-^(x,^, i) and E-^ {x,£,,t), which are 
defined as 



+E'>ix, ^, t) u{p', s') ^^5^ «(P, 5), (2) 

where [zi,Z2] denotes the gauge-link, and the renormalization scale dependence is not indicated for brevity. The 
light-like vector ri^ satisfies n{p' + p) = 2, and the skewcdness parameter ^ is defined as nA = —2^. To be specific, 
the form factors in Eq. (1) are related to the second Mellin moments of the unpolarized GPDs in (2) through [3] 

1 

dxxJ2HH^,tt)= M§ it) + ^df (t) e (3) 
-1 / 
1 

jdxxY^ Ef{x, t) = 2jQ{t) - M^{t) ~ ^ d?(i) e ■ (4) 
-1 } 

Adding up Eqs. (3) and (4) one recovers the spin sum rule [3] promising to access J'^(O), i.e. the total (spin-|-orbital 
angular momentum) contribution of quarks to the nucleon spin, through the extraction of GPDs from hard exclusive 
processes and extrapolation to the unphysical point t = 0. The sensitivity of different observables to the total angular 
momentum of in particular the u-flavour were exposed in the model studies [16, 57]. For gluons there are definitions 
and expressions analog to (2, 3, 4). Eqs. (3, 4) are special cases of the so-called polynomiality property of GPDs [14] 
stating that the A^th Mellin moments of GPDs arc polynomials in even powers of ^ of degree less or equal to A^: 



1 



h^N^iit) e^-^ for AT odd. 



/d..--i.(.,c,.)=er^(t)+4"^we^+...+J (6) 

J I Sn-iW^ for A^ odd, 

where fiavour indices are suppressed for brevity. For a spin i particle the coefficients in front of the highest power in 
^ for even moments A^ are related to each other and arise from the so-called £)-term D'^{z,t) with z = x/£, [58, 59], 
which has finite support only for \x\ < \^\, according to 

1 

/i^W(t) = -e^^(^)(i)= Uzz''~'D<^{z,t) . (7) 
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The form factors of the EMT in Eq. (1) can be interpreted [19] in analogy to the electromagnetic form factors [21] 
in the Breit frame characterized by A*^ = 0. In this frame one can define the static energy- momentum tensor for 
quarks (and analogously for gluons) 

T^.ir, = ^ / ^ exp(zAr) {p' , S'\fQ{0)h S) (8) 

with the initial and final polarization vectors of the nucleon S and S' defined such that they are equal to (0, s) in the 
respective rest-frame, where the unit vector s denotes the quantization axis for the nucleon spin. 

The components of T^.(r, s) and e''''^rjT^(r, s) correspond respectively to the distribution of quark momentum 
and quark angular momentum inside the nucleon. The components of (T^ — ^dikT^){r,s) characterize the spatial 
distribution of "shear forces" experienced by quarks inside the nucleon. The respective form factors are related to 
T;?.(r,s)by 

jQ(i) + |jQ'(t) = jd^ve-^e^^^s^r.Tgiv^s), (9) 

rf?W + f + = /d^re--T^(r) (^V-^^-) , (10) 

M2{t) - ^ (^A/2(0 - 2J{t) + ^di(t)) = JI^ j d're— ^Too(r,s) , (11) 

where the prime denotes derivative with respect to the Mandelstam variable t. Note that for a spin-1/2 particle only 
the T°^-components are sensitive to the polarization vector. Note also that Eq. (11) holds for the sum Too = Tqq-\-Tqq 
with M2{t) = M^it) + M2{t) and J{t) and di(t) defined analogously, but not for the separate quark and gluon 
contributions - since otherwise the form factor c{t) would not cancel out. 

The form factor M2 (i) at i = can be connected to the fractions of the nucleon momentum carried respectively by 
quarks and gluons. This can be seen most conveniently by considering (1) in the infinite momentum frame, and one 
obtains 

1 1 
M§{Q) = Jdx xifl + /f)(x) , (0) = jdx xfiix) , (12) 

"J 

where fi{x) ~ H'^{x,0,0) are the unpolarized parton distributions accessible in inclusive deeply inelastic scattering. 

The form factors M^''^(<), J'^''^(i) and df''~^{t) are renormalization scale dependent (the indication of the renor- 
malization scale fi is suppressed for brevity). Their quark+giuon sums, however, are scale independent form factors, 
which at t = satisfy the constraints. 



M2{0) = ^ yd3rToo(r,s)-l, 
J(0) = Jd^re'^''s,r,Tok{r,s) = ^, 



'^i(O) = / d^r r„ (r) ( rV^ - V ) = (13) 



which mean that in the rest frame the total energy of the nucleon is equal to its mass, and that the spin of the 
nucleon is 1/2. The value of di is not known a priori and must be determined experimentally. However, being a 
conserved quantity it is to be considered on the same footing as other basic nucleon properties like mass, anomalous 
magnetic moment, etc. Remarkably, di determines the behaviour of the D-term (and thus the unpolarized GPDs) in 
the asymptotic limit of renormalization scale /i — > 00 [16]. 

The form factor di{t) is connected to the distribution of pressure and shear forces experienced by the partons in 
the nucleon [19] which becomes apparent by recalling that Tij(r) is the static stress tensor which (for spin and 1/2 
particles) can be decomposed as 

r,,(r) = s{r) - +p{r)6., . (14) 
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The functions p{r) and s(r) are related to each other due to the conservation of the total energy-momentum tensor 
by the differential equation 

28p^2s^^dj^^^ (15) 
6 or r or 

Hereby p{r) describes the radial distribution of the "pressure" inside the hadron, while s(r) is related to the distribution 
of the "shear forces" [19]. Another important property which can be directly derived from the conservation of the 
EMT is the so-called stabihty condition. Integrating Jd^r r'^(ViT*-') = by parts one finds that the pressure p(r) 
must satisfy the relation 

dr r'^p{r) = . (16) 

Further worthwhile noticing properties which follow from the conservation of the EMT are discussed in App. B. Here 
we only mention that one can express c?i(t) in terms of p{r) and s{r) as (notice the misprint in Eq. (18) of [20]) 

d.it) = 5A/./d3rM!:^,(,) = i5M.|d^rM!V:^^(,), 

dx = -\^In j d^r s{r) = ^ Mn j d^r r^ p{r) . (17) 

Let us review briefly what is known about di. For the pion di^-n can be calculated exactly using soft pion theorems, 
and one obtains c?^^ = —M^^ [58]. For the nucleon the large- A^^ limit predicts [16] 

\dl + d'i\ = 0{Nl) » \dl^di\^0{N,) (18) 

which is in agreement with lattice QCD [51-53]. The constant = + df is found negative on the lattice [51-53]. 

From model calculations in the CQSM it was estimated that d^ « —4.0 at scales of few GeV^ [38, 60]. n a simple 
"liquid drop" model di is related to the surface tension of the "liquid" and comes out negative [19]. Such a model is 
in particular applicable to large nuclei. Predictions for the behaviour of the cross section of deeply virtual Compton 
scattering off nuclei made on the basis of this model [19], have been confirmed in practical calculations assuming 
realistic models for nuclei [61]. In particular, also the D-terms of nuclei were found negative [61]. 

Finally, let us discuss an interesting connection of the constant di and the mean square radius (r|,) of the trace of 
the total EMT operator. Due to the trace anomaly [62-65] the latter is given by 

= ^ F^'^F^, + (1 + 7™) ™a^a^a ■ (19) 

For notational simplicity let us introduce the scalar form factor F{t) 

(p'|f/(0)|p) = Mn u{p')u{p) Fit) . (20) 
which can be expressed in terms of the form factors in (1) as 

F{t) = AMt) + ^ (^2J{t) - AUt)^ ^ d,{t) . (21) 

It satisfies F(0) = 1 and its derivative at t = defines the mean square radius of the EMT trace operator 

(r|)=6F'(0) = 6(^M^(0)-|^) . (22) 

Analogously, one may define the mean square radius (r|;) of the energy density operator T^^ for which one finds from 
Eq. (11) the following result 

<'^^>=<^^^(«)-5^)- 
Exploring (23) we see that (r|,) is related to the mean square radius of the energy density (r|;) as follows 



Since di is observed to be negative, one has (r|,) > {r\). 



(4) = {rl) - — i- . (24) 
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III. THE NUCLEON AS A CHIRAL SOLITON 



The effective theory underlying the CQSM was derived from the instanton model of the QCD vacuum [66, 67] which 
assumes that the basic properties of the QCD vacuum are dominated by a strongly interacting medium of instantons 
and anti-instantons. This medium is diluted with a density proportional to (pav/-Rav)* where Pav ~ 0.3 fm is the 
average instanton size and i?av the average instanton separation. It is found Pav/^av ~ [66-68], see [69] for reviews. 

Due to interactions with instantons in this medium light quarks acquire a dynamical ( "constituent" ) quark mass 
which is strictly speaking momentum-dependent, i. e. M = M{p), and drops to zero for momenta p 3> Pav^. At low 
momenta below a scale set by « 600 MeV the dynamics of these effective quark degrees of freedom is governed 
by the partition hmction [70, 71] 

Vi^V'i^VU el^■p{iS<,f^{iJ,1l>,U)^ , Ses{iJ,ip,U) = J d^x {i^ - M W"^ - m)il' . (25) 

Here we restrict ourselves to two light flavours, U = exp(ir''7r") denotes the chiral pion field with U"*'^ = exp(i75T°7r''), 
and m = m,i = nid is the current quark mass neglecting isospin breaking effects. The smallness of the instanton packing 
fraction Pav/Rav plays an important role in the derivation of (25) from the instanton vacuum model. 

In practical calculations it is convenient to replace M{p) by a constant mass M = A/(0) — 350 MeV following 
from the instanton vacuum [69], and to regularize the effective theory by means of an explicit (e. g. proper-time, or 
Pauli-Villars) regularization with a cutoff of 0{p~^) whose precise value is fixed to reproduce the physical value of 
the pion decay constant /^r = 93 MeV given by a logarithmically UV-divergent expression in the effective theory (25). 
For most quantities the effects of different regularizations are of 0{M^p^^) oc p%j/R%j, i.e. parametrically small. 

The CQSM is an application of the effective theory (25) to the description of baryons [25, 26]. The Gaussian 
integral over fcrmion fields in (25) can be solved exactly. The path integral over pion field configurations, however, 
can be solved only by means of the saddle-point approximation (in the Euclidean formulation of the theory). This 
step is strictly justified in the large-A^c limit. In the leading order of the large-A^c limit the pion field is static, and 
one can determine the spectrum of the one-particle Hamiltonian of the effective theory (25) 



H\n) 



H 



(26) 



The spectrum consists of an upper and a lower Dirac continuum, distorted by the pion field as compared to continua 
of the free Dirac-Hamiltonian Hq (which follows from H in (26) by replacing U'^'^ 1) and of a discrete bound state 
level of energy E'lev, if the pion field is strong enough. By occupying the discrete level and the states of the lower 
continuum each by Nc quarks in an anti-symmetric colour state, one obtains a state with unity baryon number. The 
soliton energy Egoi is a functional of the pion field 



lev 



{En — E„ 



E„<0 



(27) 



^'soi[C^] is logarithmically divergent, see App. X for the explicit expression in the proper-time regularization. By 
minimizing Esoi[U] one obtains the self-consistent solitonic pion field Uc. This procedure is performed for symmetry 
reasons in the so-called hedgehog ansatz 



tt" (x) = e"^ P{r) , [/(x) = cos P{r) + ir^e", sin P{r) 



(28) 



with the radial (soliton profile) function P{r) and r = |x|, e- 
self-consistent profile satisfies Pc(0) = — tt and behaves as 

A 



x/r. The nucleon mass A/jv is given by Eso\\Uc\- The 



P{r) 



(1 + niT^r) cxp(— mTT?") at large r, with A 



3.9A 



(29) 



where gA = 1-26 is the axial coupling constant and the pion mass is connected to m in (25) by the Gell-Mann- 
Oakes-Renner relation for small to. In the large- A^c limit the path integral over U in Eq. (25) is solved by evaluating 
the expression at Uc and integrating over translational and rotational zero modes of the soliton solution in the path 
integral. In order to include corrections in the 1 /A^c-expansion one considers time dependent pion field fiuctuations 
around the solitonic solution. In practice hereby one restricts oneself to time dependent rotations of the soliton field 
in spin- and flavour-space which arc slow due to the large moment of inertia of the soliton, / = 0{Nc), given by 



6 ^ 



(n|T"|m) (TO|r"|n) 



En 



E„ 



(30) 



rcg 



As indicated, / is logarithmically divergent and has to be regularized. In (30) the sum goes over occupied ("occ" 
states n which satisfy £"„ < -Eicv, and over non-occupied ("non") states to which satisfy Em > S'lcv 
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IV. FORM FACTORS OF THE ENERGY MOMENTUM TENSOR IN THE CQSM 



The gluon part of the EMT is zero in the effective theory (25), since there are no exphcit gluon degrees of freedom. 
Consequently in the model the quark energy momentum tensor is conserved by itself, and the form-factor c{t) in 
Eq. (1) vanishes. This is demonstrated explicitly in App. C. The nucleon matrix elements of the effective operator 
for the quark energy momentum tensor (we omit in the following the index Q) is given by the path integral 

{p'\f^AO)\p) = ^lim^^^yd^xdV e*P'y-'P-yPV'2?V52?t/^iV'(-r/2,y)f;f(0)ji,(T/2,x)exp(^*5eff(V;,V',f/)) , 

(31) 

where JNix) denotes the nucleon current, see [26, 28, 29] for explicit expressions. The symmetric energy momentum 
tensor for quarks in the effective theory (25) is given by (the arrows indicate on which fields the derivatives act) 



(32) 



For the calculation of the EMT nucleon matrix elements in the model we have to evaluate consistently the nucleon- 
bispinor expressions appearing on the right-hand-side (RHS) of Eq. (1) in the large- A'^^ limit where = Afjy = 0{Nc) 
and Ip'I = 0{N^) such that |i| ^ Mf^. Keeping in mind that in the large-A^c limit the form-factors behave as [16] 



M2{t) = OiN^,) , Jit) = 0{N^) , d,it) = 0{N^) , 
we obtain from (1) the following relations for the form factors 

{p',S':,\t:^\p,Ss) = 2ds's(A^A''-S'''AAd,{t) 



{p',S',\T^i\p,S,) 



(33) 

(34) 
(35) 
(36) 



The expressions for M2{t) and di{t) could, of course, be separated which we shall do more conveniently at a later 
stage. Evaluating the respective components of the EMT in (31) yields (vacuum subtraction analog to (27) is implied) 



{p', 5^|f;f |p, 53) = 2M^iV, I ('^l 7°7'p'=} + (» ^ k)) \n) 



Em. — 



(37) 
(38) 
(39) 



These expressions are logarithmically divergent and have to be regularized appropriately, see App. D for details. 

Noteworthy, the matrix elements for the components T^^ and T^^ related to M2{t) and di{t) are spin-independent 
and receive contributions from leading order of the large-A'c expansion. In contrast to this, in order to address T^^ 
connected to J{t) one needs matrix elements involving nucleon spin flip which appear only when considering 1/A^c 
("rotational") corrections. Inserting the results (37, 38, 39) into Eqs. (34, 35, 36) yields 



M2{t) 



diit) 
di{t) 
J{t) 



M, 



N 



J2 i?„(n|e*^*|n) 



n.occ 



n,occ 



rcg 



8It 



E 



{n\r'\j) 



Ej - En 



(40) 
(41) 
(42) 



rcg 
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The derivation of Eqs. (40, 41, 42) follows closely the derivation of the model expressions for electromagnetic [28] or 
other form factors, and we omit the details here. Instead, we demonstrate explicitly in the Appendices E-H that one 
obtains the same expressions for the form factors from the model expressions for GPDs via the sum rules (3) and (4). 
We introduce the Fourier transforms of the form factors which are radial functions ("densities") defined as 



p{r) 



Nc En rnir)Mr) 



reg 



Nr. 



J2ct>*n{r) (7°7P)0n(r) 



reg 



n.occ 



24/ 



(43) 
(44) 
(45) 



where it is understood that jp = |(^"' — ^"'), and which allow to reexpress M2{t), J{t) and di{t) in (40, 41, 42) as 



M2{t) 



di{t) = ^ J d'r pEir) joirV~t) 



diit) 



15M. 



d r p{r) 



J{t) = 3 / d^r pj{r) 



3i{rV-t) 



(46) 
(47) 
(48) 



Here jk{z) denote Bessel functions with jo{z) — and ji{z) — — jQ(z). 

The densities (43, 44, 45) are convenient not only because their numerical evaluation is more economic than the 
direct calculation of the form factors (40, 41, 42). These densities are interesting objects by themselves, and it is 
instructive to discuss their theoretical properties in detail which we shall do in the following. Simultaneously we will 
present the numerical results for the densities. 

For the numerical calculation we employ the so-called Kahana-Ripka method [72] , whose application to calculations 
in the CQSM is described in detail e.g. in Ref. [29], and use the proper-time regularization. The latter allows to include 
effects of symmetry breaking due to an explicit chiral symmetry breaking current quark mass m in the effective action 
(25), and to study in the model how observables vary in the chiral limit. (For a study — in the spirit of [55] — of 
observables at pion masses as large as they appear in present day lattice calculations the reader is referred to [56].) 

In order to explore effects of different regularizations we perform also a calculation with the Pauli-Villars regular- 
ization method which, however, is applicable only in the chiral limit [73]. All results are summarized in Table I. 



ENERGY DENSITY 



The density pEir) is just the energy density T'^'^(r) in the static energy momentum tensor (8). By using the 
ortho-normality of the single-quark states Jd'^r 0* (r)0„/(r) = Snn' and comparing to Eq. (27) we find that 



j d'^v pE{r) ^ N, Y.En j 



d3r0;(r)0„(r) 



N 



(49) 



The normalization (49) ensures the correct constraint of the form factor M^it) at t = 0. In fact, by taking in Eq. (46) 
the limit i ^ (notice that di{t) takes a well-defined finite value for t ^ 0, see below) we obtain 



M2(0) = ^ yd'r pE{r 



= 1 



(50) 



This is the consistent constraint in the model for Ahit) at t = 0, cf. Eq. (12), since there are no gluons in the effective 
theory such that consequently the entire momentum of the nucleon is carried by quarks and antiquarks [32] . 

Fig. la shows the normalized density pE{r)/M]y as function of r for the physical situation with a pion mass of 
140 MeV. In this case in the model the nucleon mass is about 1250 MeV. This overestimate of the physical nucleon 
mass of 0(300 MeV) is typical for the soliton approach and its origin is well-understood [74]. In the center of the 
nucleon one finds p_e(0) = 1.7GeVfm~^ or 3.0 x 10^^ g cm~^. In order to gain some intuition about this number we 
remark that this corresponds roughly to 13 times the equilibrium density of nuclear matter. 
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FIG. 1: (a,) The normalized energy density pE{r) / Miq of tlie nucleon as function of r for the physical situation with m,r = 140 MeV. The 
curve is normalized such that it yields unity upon integration over the entire volume, (b) The normalized energy density Avr^ pE{r)/Mi\f 
as function ofr in the chiral limit and for m,r = 140 MeV. The curves arc normalized such that one obtains unity upon integration over r. 



It is instructive to consider the energy density in the chiral hmit. The resuh is shown in Fig. lb where we compare 
47rr^P£(r)/Mjv as functions of r for TOtt = and 140 MeV. The curves are normahzcd such that one obtains unity 
when integrating over r. Fig. lb shows that with decreasing TTT-Tj- tliG GiiGrgy d-Giisity is sprecicl morG widely. This Cctn 
be quantified by considering the 77171-- dependence of the mean square radius of the energy density (23) defined as 

. . (51) 

J d'^r pE[r) 

and which increases in the chiral limit, see Table I. This is an intuitively expected feature. As the pion mass decreases, 
the range of the "pion cloud" increases and the nucleon becomes "larger" . 

The popular idea of the nucleon consisting of a "quark core" surrounded by a "pion cloud" is strictly speaking 
well defined in models only. Here, in the CQSM, we shall associate the contribution of the discrete level as "quark 
core" and the contribution of the negative continuum states as "pion cloud" . From the long-distance behaviour of 
the soliton profile (29) one finds in the chiral limit 



2 



PEir) = ^ trp 



V'=t/(r) V'=C/t(r) 



largo / 3gA \ 1 



„6 



(52) 



where tr^ is the trace over flavour indices of the SU(2) matrices and the dots in the intermediate step denote terms 
which contain higher U -field derivative terms and vanish faster at large r than the displayed leading term. The result 
(52) can be read off from Eq. (7.8) of [32]. At ttItt ^ the decay of psir) at large r is exponential due to the 
corresponding behaviour of the soliton profile (29). This diminishes the "range of the pion cloud" and reduces {r'^)- 
These observations can be further quantified by considering the chiral expansion of (51) which gives, see App. I, 

/.2^_/02^ 81gl _ ^ ^^^^ 



647r/2Mw 



Here and in the following the ° above a quantity denotes its value in the chiral limit, and the dots denote terms 
vanishing faster in the chiral limit than the respective leading term. Considering the non-commutativity of the limits 
Nc ^ oo and (see the discussion below in Sec. VIII) Eq. (53) agrees with chiral perturbation theory [45]. 

The term linear in in Eq. (53), i.e. the leading non-analytic (in the current quark mass m cx m^) contribution 
to the mean square radius of the energy density accounts almost entirely for the reduction of (r|;) from TOtt = to 
140 MeV, see Table I.^ At the physical pion mass we find {r%) = 0.67 fm^. This value is similar to the electric charge 
radius of the proton. In fact, we observe a qualitative similarity of the energy density and electric proton charge 
distributions in the model [28, 29]. 



^ In the CQSM the physical value of qa is underestimated by about 30% in the proper time rcgularization in the leading order of the 
1/Nc expansion to which we work here. For consistency we have to use here this leading order Nc model result for gA- Including 1/Nc 
corrections the model describes qa more accurately [29]. 
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FIG. 2: (a.) The normalized angular momentum density pj(r)/Ji^ as function of r for the physical situation with m,r = 140 McV. The 
curve is normalized such that it yields unity upon integration over the entire volume, (b) The normalized angular momentum density 
4:nr^ p J {r) / J ff as function of r for m-n- = and 140 McV. The curves arc normalized such that one obtains unity upon integration over r. 



VI. ANGULAR MOMENTUM DENSITY 



Taking in Eq. (48) the limit t yields 

J(0) = / d^r pj{r) , (54) 



which shows in which sense it is adequate to refer to pj{r) as the "angular momentum density". In order to see that 
J(0) = ^, i.e. that the constraint (13) is satisfied in the CQSM we insert (45) in the above equation and obtain 

„,oco J n reg _^ f 

(55) 

In the intermediate step in Eq. (55) we recovered the model expression for the second moment of X]/ (H^ + E^){x, ^, t) 
at i = 0, see Appendix C of [42]. This sum rule, which follows from adding up Eqs. (3) and (4), was explicitly proven 
to be satisfied in the model in [42]. In the model the entire nucleon spin is due to the spin and orbital angular 
momentum of quarks and antiquarks. and hence 2J(0) = 1 [42]. This again is a correct and consistent result since 
there are no explicit gluon degrees of freedom in the CQSM. 

The numerical result for the normalized angular momentum density p.j{r)/JN = '2-Pj{i') as function of r for the 
physical situation is shown in Fig. 2a. {Jm = \ denotes the nucleon spin.) We observe that pj{r) oc at small r. 

In Fig. 2b we compare the normalized angular momentum densities Airr^ pjir) / Jjs[ as functions of r for ttItt = 
and 140 MeV. These curves are normalized such that one obtains unity upon integration over r. Within the rotating 
soliton picture the result is reasonable. The smaller tjIt^^ the larger the nucleon, and the more important is the role 
of the region of large r for the description of the soliton rotation, i.e. for the spin structure of the nucleon. This is 
reflected by the mean square radius of the angular momentum density which we define in analogy to (51) as 

. SffPf) . (56) 
Jd^rpj{r) 

The mean square radius of the angular momentum density increases with decreasing m7r, see Table I. In the chiral 
limit pj{r) oc ^ at large r, i.e. (rj) diverges in the chiral limit. As a consequence J{t) has an infinitely steep slope 
at t = in the chiral limit, see Sec. VIII and Appendix I for further discussions. 



10 



p(r) inGeVfm"^ (a) s(r) in GeV fm"^ (b) p(r) & s(r) in y R (c) 




0.5 1 1.5 rinfm 0.5 1 1.5 rinfm 1 rinR^ 



FIG. 3: (a) The pressure p{r) as function of r for = and 140 MeV. (b) The same for the function s{r) defined in Eq. (14) which 
describes the shear forces in the nucleon, and is related to p{r) by the relation (15). (c) p{r) and s(r) in a liquid drop in units of -f/Rd as 
functions of r in units of R^. Here R^ is the radius of the drop, and 7 is the surface tension. The 5-functions in p{r) and s{r) in Eq. (60) 
are smeared for better visibility, see text. 



VII. PRESSURE, SHEAR FORCES, SOLITON STABILITY, AND SIGN OF di 

For the pressure (44) the analogon of the "normahzation relations", Eqs. (49) and (55), of the other densities is the 
stabiUty criterion (16). Integrating the r^-weighted model expression for p(r) over r wc obtain 

00 

fdrrMr)^{^ /d^x J] ^^(x) (7%?) 0„(x) = ^ J] (n| 7°7P |n) = , (57) 

Q n, occ n, occ 

because more generally the tensor K^^ defined as 

n, occ 

is zero — however, if and only if, one evaluates the expression (58) with the self-consistent profile, i.e. with that profile 
which minimizes the soliton energy (27). This was proven in Ref. [32]. Thus, the stability criterion (16) is satisfied. 

Next let us check that one obtains in the model the correct result for the constant di = di{0) defined in Eq. (13). 
Expanding the expression (47) around small t we obtain 

d,(i) = iM^^|d3rp(r)Q + ^ + 0(t)). (59) 

Due to (57) the 1/t-term drops out, and we verify in the CQSM the relation (17) between di and the pressure. 

Fig. 3a shows the pressure p{r) as function of r. In the physical situation p(r) takes its global maximum at r = 
withp(O) = 0.23GeV/fm^ = 3.7-1034 Pa. This is 0(10-100) higher than the pressure inside a neutron star [75]. Then 
p{r) decreases monotonically — becoming zero at tq = 0.57 fm — till reaching its global minimum at min = 0.72 fm, 
after which it increases monotonically remaining, however, always negative. The positive sign of the pressure for r < tq 
corresponds to repulsion, while the negative sign in the region r > tq means attraction. 

In Fig. 3a we see how the pressure depends on the pion mass. The pressure in the center of the nucleon increases as 
increases — obvious consequence of the fact that the (energy) density also increases, see Table I. As a response 
to the increased p{r) at small r — keep in mind the stability condition (16) — the pressure takes also larger absolute 
values in the region r > rp where it is negative. This can again be intuitively understood because a heavier particle 
is more tightly bound, i.e. the attractive forces are stronger. The zero of p{r) moves towards smaller values of r with 
increasing TO^r, see Table I. 

In Fig. 3b wc show the distribution of the shear forces s(r) obtained from our results for p{r) by solving the 
differential equation (15). The distribution of shear forces is always positive. It reaches for 111^^ ~ 140 MeV a global 
maximum at r = 0.40 fm. The position of the maximum is weakly dependent on m^r. At small r we observe s{r) cx . 

It is interesting to compare to which extent the nucleon "resembles" a liquid drop of radius Rd with constant density 
and constant pressure pq inside. In such a drop the pressure and shear forces are given by [19] 

p(r) = po e{Rd ~r)-^ poRd S{Rd ~ r) and s(r) = 7 5{Rd - r) , (60) 
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r^p(r) in GeV fm"^ (a) r''p(r) x ^M^ 47i in fm"^ (b) 




0.5 1 rinfm 0.5 1 1.5 2 2.5 3 3.5 4 rinfm 



FIG. 4: (a) r^p(r) as function of r from the CQSM at the physical value of . The shaded regions have — within the numerical accuracy 
of about half percent — the same surface areas. This shows how the stability condition dr r^p(r) = in Eq. (57) is realized, 
(b) The same as (a) but with an additional power of and the prcfactor 5nMff. Integrating this curve over r yields di according to (17). 
The plot shows that one obtains a negative sign for di as a consequence of the stability condition (16) shown in Fig. 4a. 



where 7=5 PoRd denotes the surface tension. We show this situation in Fig. 3c — where, however, for better 
visibiHty the (5-functions in (60) are smeared out. This corresponds to allowing the density in the drop to decrease 
continuously from its constant inner value to zero over a finite "skin" (of the size ^ Jo-^d in Fig- 3c). 

Comparing the liquid drop picture to the results from the CQSM we observe a remote qualitative similarity. In 
contrast to the liquid drop, the density "inside" the nuclcon is far from being constant, see Fig. la, and one cannot 
expect the pressure in the nucleon to exhibit a constant plateau as in the liquid drop. Still the pressure exhibits the 
same qualitative features. The shear forces become maximal in the vicinity of what can be considered as the "edge" 
of the object. This is the case in particular for the liquid drop However, the "edge" of nucleon is far more diffuse, 
and the distribution of shear forces s{r) is widespread. Of course, the nucleon can hardly be considered a liquid drop. 
Such an analogy might be more appropriate for nuclei [19]. Nevertheless this comparison gives some intuition on the 
model results — in particular, about the qualitative shape of the distributions of pressure and shear forces. 

Next let us discuss how the stability condition (57) is satisfied. Fig. 4a shows r'^p{r) as function of r. The shaded 
regions have the same surface areas but opposite sign and cancel each other — within numerical accuracy 



dr r^p{r) ^ 2.61 MeV, 



d?' r'^p{r) 



-2.63 MeV. 



(61) 



In order to better understand how the soliton acquires sta- 
bility, it is instructive to look in detail how the total pressure 
is decomposed of the separate contributions of the discrete 
level and the continuum contribution. Fig. 5 shows that the 
contribution of the discrete level is always positive. This con- 
tribution corresponds in model language to the contribution 
of the "quark core" and one expects a positive contribution 
("repulsion") due to the Pauli principle. At large r the dis- 
crete level contribution vanishes exponentially since the dis- 
crete level wave- function does so [26] . 

The continuum contribution is throughout negative — as 
can be seen from Fig. 5 and can be understood as follows. 
The continuum contribution can be interpreted as the effect 
of the pion cloud which in the model is responsible for the 
forces binding the quarks to form the nucleon. I.e. it pro- 
vides a negative contribution to the pressure corresponding 
to attraction. In the chiral limit the continuum contribution 
exhibits a power-like decay which dictates the long-distance 
behaviour of the total result for the pressure as follows 



r^p(r) in GeV fm'^ 



0.05 - 



-0.05 - 



quark core 

pion cloud 

total 



0.5 



1,5 r in fm 



FIG. 5: Tile pressure p(r) as function of r for m,r = 140 MeV. 
Dotted line: Contribution of the discrete level associated with 
the quark core. Dashed line: Continuum contribution associ- 
ated with the pion cloud. Solid line: The total result. 



p{r) = - 



3gA 



— and s{r) =3 



— at large r. 



(62) 
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where for completeness we quote also the result for s(r). For ^ the continuum contribution exhibits an 
exponential decay at large r due to the Yukawa tail of the soliton profile (29) — like the contribution of the discrete 
level, however, it still dominates the behaviour of the total result. Fig. 5 reveals that the actual cancellation between 
the different contributions leading to (57) are even more impressive than we could guess considering the numbers in 
Eq. (61). For = 140 MeV we obtain the following numbers, and see that the two contributions cancel each other 
within numerical accuracy 



oo oo 

/ dr r^p(r) = 26.69 MeV , / dr r^p(r) 

lev ^ 



-26.71 MeV. (63) 



Finally, we discuss the relation of stability and the sign of the constant di. By comparing Figs. 4a and 4b we 
immediately understand that in the CQSM the constant di takes a negative value 

di < . (64) 

This result is presumably of general character. In fact, the observation (64) follows naturally from our intuition on the 
pressure and shear forces distributions we gained from our study. It seems physically intuitive that in a mechanically 
stable object the following conditions hold 

{> for r < To 
= for r ~ ro 
< for r > ro 



(II) distribution of shear forces satisfies for all r s(r) > . 

From (I) and the stability criterion (16) we conclude that di cx J^dr r'^p{r) < as illustrated in Figs. 4a and 4b. 
The same conclusion follows from (II) using di oc ~ J^dr r^s{r) < in (17). 

Of course, one can imagine a pressure distribution with more zeros and a different shape than in (I) still yielding 
(64). However, (I) is the simplest case which one may expect to hold for a ground state — like the nucleon. Such a 
ground state object is characterized by having one "surface" only, although a quite smeared out one in the case of 
the nucleon, hence the condition (II). A more complicated distribution of the pressure — with more zeros, i.e. also 
more maxima and minima — would imply an object with several (smeared out) surfaces, which follows from the last 
condition in App. B. 

The conjecture (64) is — besides being physically appealing — in agreement with all available information on di, 
see Sec. II. One may therefore suspect that (64) is a general theorem which connects the stability of an object to the 
sign of its constant di. However, such a theorem — if it exists — remains to be rigorously proven for the general case. 
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VIII. RESULTS FOR THE FORM FACTORS 



From the densities p£;(r), p.j{r) and p{r) which we discussed in detail in Sees. V, VI and VII we obtain the form 
factors of the EMT by means of Eqs. (46, 47, 48). The results are shown in Fig. 7 for = and 140 MeV. In the 
CQSM the form factors M2(t), J{t) are normalized at t = as M2(0) = 2J(0) = 1 as proven in Sees. V and VI. The 
numerical results satisfy these constraints within a numerical accuracy of (1-2)%, see Figs. 7a and 7b. 

In contrast, the normalization of the form factor di[t) at t = is not known a priori. We find di ~ di{0) < as 
anticipated in Sec. VII. The constant di has a well-defined chiral limit. This can be concluded from Eq. (17) which 
relates di to the distributions of pressure or shear forces, and from the fact that p{r) and s(r) drop off sufficiently 
fast at large r in the chiral limit, see (62). 

The absolute value of di decreases with increasing and 
we observe a strong sensitivity of di to m^. That this is 
not surprizing can be understood by considering the chiral 
expansion of di. Expanding di(m.,r) in the model for small 
m,r we obtain, cf. Appendix I, 







dl(mjr) 



-1 - 



(ii(m^) ^ di + 



64 7r/2 



(65) 



Thus, we see that di (m,r) receives a large leading non-analytic 
contribution in the current quark mass m cx , which is the 
origin of the strong m^-dependence of diirriTr)- The expan- 
sion (65) approximates in the model di(m7r) to within an 
accuracy of 15% up to the physical point. At the physical 
point the absolute value of di is reduced by about 30 % with 
respect to its chiral limit value. The numerical results for di 
are summarized in Table I. 

Next, let us focus on the derivatives of the form factors at 
t — Q. From Eq. (11) we conclude that at t = the slope of 
the form factor M2{t) is given by 



exact model result 

chiral expansion up to 0(m^} 



_1_ 



50 



100 140 



nijj in MeV 



FIG. 6: The constant d\ as function of m,r. Comparison of 
the full CQSM resuJt, and the expansion of di{m^) up the 
leading non-analytic term according to Eq. (65). The square 
marks the physical point. 



Mm 



{r 



6 



di 



A/ 2(0) 



21 

128^/2 



+ 



N 



(66) 



where the chiral expansion follows from (53) and (65). In particular, ^^2{^) is well-defined for all including the 
chiral limit. The situation is different for the form factors J{t) and di{t). The slope of J{t) at t = can be deduced 
from (48) and is given by 



J'(0) 



6 



(67) 



with (rj) defined in (56). Since the mean square radius of the angular momentum density diverges in the chiral limit, 
see Sec. VI, so does J'(0). See also the discussion in App. I. 

That also the slope of di{t) becomes infinitely steep at t = can be understood as follows. From Eqs. (11, 14, 15) 
it follows that the derivative of di{t) at < = be expressed as 



d[{t) 



t=o 



Mn f ,3 4 . . Mn f 3 4 , , 

'l2~ J "^^"^ " ^ Jdrr p{r) . 



(68) 



Since in the chiral limit p{r) and s(r) drop off as -pi at large r we see that d{ (0) diverges as 0. 

To make this statement more quantitative let us expand di{t) in the chiral limit in the model for small t. We obtain 



di (t) = di 



45 gjMjv 
512/2 



-t + 0{t) , 



(69) 



which means that in the chiral limit d{ (t) (x 1/^/^ at small t. Alternatively, one may keep ttItt 7^ 0, evaluate the 
derivative of di{t) at t = 0, and then consider small m,r- (Note that the limits 111^^ —>■ and t ^ do not commute.) 
Then we find that the slope of di{t) at t = diverges in the chiral limit as 
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FIG. 7: The form factors of the energy momentum tensor M2{t), J(t) and di(t) as functions of t for the pion masses = 0, 140 MeV. 
All form factors can be well approximated by dipole fits, however, with the exception of J{t) and di{t) in the chiral limit which exhibit 
inhnitoly steep slopes at t = 0, see text. 



The results (70) and (70) are derived in Appendix I. The numerical results for J{t) and di{t) in Figs. 7b and 7c 
indicate the infinitely steep slopes at t = within the numerical accuracy. 

Notice that we derived the analytical results in (68, 69, 70) in the framework of the CQSM. However, the leading 
non-analytic terms (i.e. terms oc m.^-) in the chiral expansion of di(t) in (68, 69, 70) are dictated by chiral symmetry 
breaking, and are independent of the details of the chiral theory chosen to derive them. In fact, our results (68, 69, 70) 
agree with those obtained from chiral perturbation theory in Ref. [45] — provided one takes into account an important 
difference. The CQSM is formulated in the large- A^c limit which does not commute with the chiral limit [47]. At 
large Nc the masses of the nucleon and A-resonance are degenerated: Ma — Afjv = C>{N~^). Therefore in the CQSM 
— in addition to the nucleon considered in chiral perturbation theory [45] — the A-resonance contributes on equal 
footing as intermediate state in chiral loops. Considering that in the large A^c-limit the pion-Delta-nucleon and pion- 
nucleon- nucleon couplings are related as g-j^NA = ^QttNN (phenomenologically satisfied to a very good approximation), 
one finds that the A-resonance makes a contribution to leading non-analytic terms which is — for scalar-isoscalar 
quantities — two times larger than that of the nucleon [48]. Hence, our leading non-analytic terms in Eqs. (68, 69, 70) 
are 3 times larger than those obtained from chiral perturbation theory in Ref. [45] where Nc was kept finite.^ Other 
examples of the derivation of leading non-analytic terms in chiral soliton models can be found in [31, 48]. 

Next, we turn to the discussion of the t-dependence of the 
form factors, and recall that in the large- A'c limit we strictly F(t)/F(0) 
speaking arc restricted to \t\ ^ M^. However, in practice 
it is observed that the CQSM provides reliable results for 
electromagnetic form factors up to \t\ < 1 GeV^ [28, 29]. 

Fig. 7 shows the form factors of the EMT as functions of 
t for |t| < 1 GeV^ for = and 140 MeV. For 7^ 
the form factors can be well approximated by dipole fits 



m 



(71) 



with the values for the dipole masses quoted in Table I. 

It is instructive to compare these results within the model 
to the electromagnetic form factors — for definiteness we 
choose the electric form factor of the proton GE{t) com- 
puted in the CQSM in Ref. [28]. Fig. 8 shows that J{t) 
and Gsit) exhibit a similar t-dependence. However, M2{t) 
falls off with increasing \t\ slower than GE{t), while di{t) 
shows a faster fall off. 




t in GeV 



FIG. 8: EMT form factors as functions of t. Dashcd-dottcd 
line: M2{t). Long-dashed line: J{t). Dashed line: di{t); all 
computed here. Solid line: electric proton form factor Gsit) 
from Ref. [28]. All form factors refer to the physical pion mass, 
and are normalized with respect to their values at t = 0. 



^ Actually, there is one more subtlety to be considered. In the result for the form factor di{t) = ^C2{t) in Eq. (20) of [45], see App. A 
for the discussion of the notation, in addition the constant a^'^ = M^'^ appears, which describes the fraction of the pion momentum 
carried by quarks. In the effective theory (25) quarks carry the entire momentum of the pion, i.e. M^'"" = 1 [58]. 
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TABLE I: The pion mass dependence of different quantities computed in the CQSM: the energy density in the center of the nuclcon 
Pe{0), the mean square radii {r'^) and (r^) as defined in Eqs. (51, 56), the pressure p(0) in the center of the nucleon, the zero of the 
pressure defined as p(ro) = 0, the constant di, the dipole masses of the form factors M2(t), J(t) and rfi(t) as defined in Eq. (71), and the 
mean squared radius of the trace of the EMT (22). In the chiral limit J(t) and d\{t) have infinitely steep slopes at t = 0. In these cases 
dipole fits do not provide useful approximations and arc undefined (labelled by " — " in the Table). The proper time regularization can be 
applied, in principle, to any [55]. The Pauli-Villars regularization is applicable only for = [73]. 
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Pauli-Villars regularization: 




















0.75 


0.86 


oo 


0.332 


0.63 


-4.75 


0.804 






1.24 



One popular assumption in literature in the eontext of modelling GPDs is to assume a generic factorizcd Ansatz 
of the type H(^x,£,,t) = H{x,S,)G{t) where G{t) denotes the respective form factor (other approaches are discussed 
in [16, 76-80]). This assumption implies that the form factors of the EMT should have approximately the same 
t-dependence as the electromagnetic form factors. Our results indicate that this is quite a rough approximation, and 
support the observations that in the CQSM H{x,£^,t) ^ H{x,£^)G{t), see Refs. [16, 38] 

Let us compare our result for di with the result from direct calculations of GPDs in the model [38] which yield'^ 
-8.0 at the low model scale in the chiral limit. The discrepancy with the corresponding results in Table I is 
due to the fact that in [38] the momentum-dependent constituent mass M{p) was employed vs. our proper-time or 
Pauli-Villars regularization with M = const, and that the approximative "interpolation formula" was used vs. the 
exact numerical calculation done here. 

For the mean square radius (rp) of the operator of the EMT 
trace, see Eqs. (23, 24), we obtain in the chiral limit 



7Q 



(rp) infm^ 



= (1.0- 1.2) fm^ 



(72) 



while its value at the physical point is reduced by about 0.2 fm 
compared to (72), see Table I. From the chiral expansion of 



which follows from Eqs. (53, 65, 24) and reads 



64^/2M, 



TV 



(73) 



we see that the leading non-analytic term in (73) explains the 



' ' ' ' I ' ' ' I 

exact model re.sLdt 
chiral expansion up to O(m^) 




main portion of the observed reduction of (rp) at the physical 



in MeV 



50 100 140 

point compared to its chiral limit value, see Fig. 9. "''"^^'^ ^^^^"^ ^'"^^ "P^""*"' ™- 

It is instructive to compare our result (72) for the mean square radius of the EMT trace to the mean square radius 
of the traceless part of the EMT (0.3 fm)^ estimated by means of QCD sum rules [81]. The instanton vacuum model 
provides a possible explanation why the radius of the trace of the EMT is so much larger than the traceless part. In 
chiral limit the trace part is basically the gluonic operator F^i,F^'^ which is due to one- instanton contributions and 
appears in leading order in the instanton packing fraction. The traceless part arises from instantoii-aiiti-instanton 
contributions which appear at subleading order [68]. 



From this number it was estimated that di 



predicts for the ratio /Ml 



Q 



-4.0 at experimentally relevant scales of few GeV^ [60] in the following way. The model 



1.0 and experimentally one finds M? ~ 0.5 at scales of few GeV'^. This estimate neglects strictly 



speaking the different evolution properties of and which is justified because uncertainties in the model dominate. 
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IX. CONCLUSIONS 

111 this work we presented a study of the form factors of the EMT in the large- A^c hmit in the framework of 
the CQSM. We provided numerous checks of the theoretical consistency of the model results. Among others, we 
demonstrated that the same model expressions for the form factors can be derived from the EMT and from GPDs. 

We computed the spatial density distributions, and mean square radii of the operators of different components of 
the EMT and its trace. Interesting results are that the energy density related to T^" in the center of the nucleon is 
about 1.7 GcV '^, i.e. about 13 times higher then the equilibrium density of nuclear matter. The mean square radius 
of the operator T*"^ is about 0.67 fm^. For the mean square radius of the "angular momentum distribution" related 
to the operator {k = 1, 2, 3) we find a much larger value 1.32 fm . 

We studied the spatial distribution of strong forces in the nucleon as described in terms of the distributions of 
pressure and shear forces which are defined by the spatial components T^^ (j, fc = 1, 2, 3) of the EMT. As a byproduct 
of this study we learned how the soliton acquires stability in the CQSM — namely due to subtle balance of repulsive 
forces in the center due the "quark core" and attractive forces in the outskirts of the nucleon due the "pion cloud" 
which bounds the quarks in the model. 

We observed a physically appealing connection between the criterion for the stability of a particle, and the sign 
of the constant di, i.e. the form factor di(t) at zero-momentum transfer. Our observations imply that for a stable 
particle one always has di < 0, though we cannot prove this conjecture for the general case. All available results for 
di in literature are compatible with this observation. 

We derived the leading non-analytic chiral contributions to the form factors in the large Nc limit, which agree with 
results from chiral perturbation theory [44-46] provided one takes into account that the limits Nc oo and 
do not commute [47, 48]. 

We observed that the model results for the form factors of the EMT can, for ttItt > 0, be well approximated by dipole 
fits. The different form factors have different ^-dependences. For J^{t) we obtain in the model a dipole mass similar 
to that of the electromagnetic form factors of the proton. The dipole mass of M2{t) is much larger than that, while 
that of di{t) is substantially smaller. These results are of interest for the phenomenology of hard exclusive reactions, 
and in particular for the task of extrapolating J^{t) to t = Q which is necessary to extract from data information on 
what portion of the nucleon spin is due to quarks. 

Our results yield in the chiral limit for the constant di = —(3.5 — 4.8) depending on the regularization and confirm 
sign and — within model accuracy — magnitude of previous results [38, 60]. We observe, however, also a strong 
m,r-dependence of di which is dominated by a sizeable leading non-analytic (in the current quark mass) contribution 
proportional to mT,. The latter is responsible for reducing the absolute value of di by about 30% at the physical point 
compared to its chiral limit value. 

Wc estimated the mean square radius of the EMT trace operator to be about (Ifm)^ which appears much larger 
than the corresponding mean square radius of the tracclcss part of the EMT found to be (0.3 fm)^ [81], and noticed 
that the instanton vacuum model provides a possible explanation for that. 

A study of the m^-dependence of the EMT form factors in the model at larger values of in the spirit of Ref. [55] 
and a detailed comparison of the model results to lattice QCD data [49-54] will be presented elsewhere [56] . 
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Note added After this work was completed the work [88] appeared where in particular the constant di was studied. 
For d"^'' similar results were obtained using somehow different model parameters. Interesting is the estimate for the 
flavour combination d\'^'^ which was found rather small and confirms the large- A'c prediction (18). 
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APPENDIX A: ALTERNATIVE DEFINITION OF FORM FACTORS 

The following alternative definition of form factors of the EMT is commonly used in literature, see e.g. [3], 



(p'|T;?.^(o)b) = Hp') 



P 



■ b)- (Al) 



By means of the Gordon identity 2Mmu'^°'u — u' {icr°"^ IS.^ + 2P")u Eq.(Al) can be rewritten as Eq.(l) with 

C^'^{t) = -d'^'^it). (A2) 

The constraints (13) translate in this language into A'^ (0) + A'^ (0) = 1 and BQ{0) + B'^{0) = 0. The latter constraint 
means that the total nucleon " gravitomagnetic moment" vanishes. 

In models, in which the only dynamical degrees of freedom are effective quark degrees of freedom, the constraint 
B'^{0) — must hold. Such is the situation in the CQSM where consequently this constraint is satisfied [42]. 

Interestingly, it was argued [82] that also in QCD the quark and gluon gravitomagnetic moments of the nucleon 
could vanish separately, i.e. 5*3(0) = and B'^{0) = 0. That would imply that M^{0) = 2 J'5(0) and M^{0) = 2 J<^(0) 
at any scale, and not only in the asymptotic limit of a large renormalization scale [3], see [82] for details. 

APPENDIX B: GENERAL RELATIONS FROM THE CONSERVATION OF EMT 

Here wc collect some worthwhile noticing relations for p(r) and s(r) which can be derived from the differential 
equation (15) — i.e. which follow from the conservation of the EMT. 

• For any s{r) one obtains from (15) a p{r) which automaticaiiy satisfies the stability condition (16). Therefore, 
when computing these quantities for example in a model, the computation of the pressure is more important 
and reliable in the sense that the result can be cross-checked by the stability condition (16). 

• The pressure at the origin is connected to s(r) by the following integral relation 

oo 

p{0) = 2 Jdr^ . (Bl) 



• Assume that p{r) and s(r) vanish at large r faster than any power of r to justify below integration by parts. 
(In Sec. VII we have seen that this is always the case with the exception of the chiral limit.) Then the following 
relations hold between the "Mellin moments" of s(r) and p{r) 

no oo 

drr^5(r) = -|^±l|y'drr^p(r) for iV > -1, (B2) 



which are valid also for non-integer values of N. Eq. (17) quoted in Sec. II is just a special case of (B2). 

• The case = 2 in (B2) has to be treated with care. Taking the limit ^ 2 on the right-hand-side of (B2) 
yields — upon use of (16) — the following result 

oo oo 

9 f 

dr sir) = — — dr logr p{r) . (B3) 



One may worry in which units the "r" in "log?*" should be provided in (B3). However, this is immaterial due 
to the stability condition (16). 

p(r) and r^s{r) exhibit extrema at the same r. 
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APPENDIX C: CONSERVATION OF EMT IN THE CQSM 

In this appendix we demonstrate explicitly that {N'\d^Tl^^\N) ~ 0. For that we directly evaluate in the model 
matrix elements of the operator d^T^^ . For <9oT|?ff ^ derivation analog to that yielding (37) and (39) yields immediately 

(7V'|9otffl^>=0. (CI) 

Hereby one has to consider the cases v = Q and v ^ k separately, since the respective model expressions arise from 
different orders in the 1/Nc expansion. 

A derivation analog to that yielding (38) gives for matrix elements of the operator diT^g the following result 



n.occ 



. MnN, 

OS'S. 



n,occ 



y"d3xe*^^(Ki +/\2) (C2) 



where 

K, = -C[V'^(-J7VVO0„)] + [-*7°7''VV„]*[V^(/.„] + (V'=0;)(-^7°7''VV„) - [V'^(-*7°7>„)1</>„ 

K2 - +z0:7V(V'</'«)-»(V"0;)7V0n • (C3) 

From the single quark equations of motion (26) one obtains the identities (— j7"7*V*)0„ = (En — M^^U^^)(f)n and 
(— V^0„) = [En — M"^ — iM-f^ {W^ U'''^))(f)n which allow respectively to rewrite Ki and K2 as follows 

/vi = -K2 = 2A/,/);7°(V^[/'^^)./)„ . (C4) 

Thus, {N'\dif^^\N) = in Eq. (C2). This proves that the form factor c(i) in Eq. (1) vanishes in the model. This is 
consistent since, due to the absence of gluons, in the CQSM the quark part of the EMT must be conserved by itself. 

APPENDIX D: REGULARIZATION 

The proper time regularized model expressions for the contimmm contributions to the form factors (40, 41, 42) read 



M2(t),o„t - 7T^di(t)eont = ^ E i?i (^^n , A) (n| e'^* | n) , (Dl) 



diitUnt = R2iEn,A){n\h^p, e*^4|n) , (D2) 

ri,all ^ 
~klm A 



J(i)co„t = ^-^^ Y i?3(i?„,i?„ A)(n|r'|j)(j| ( e^^^p™ \ + (K + i?,)e'^V7'" J l^i) (D3) 

TLjj all \ ^ ^ / 

with the regulator functions Ri defined as follows 



n / A \ 1 [da 2\ 
^^^""'"^^^ 47^ j ^ exp(-au. ) 

1/A2 

i?2(c^,A) = -^ |ji?i(^,A) 



CX3 



i?3(wi,tj2, A) = y"da 

1/A2 



exp(— awj) — exp(— q;cj|) wi exp(— aw^) + 0^2 exp(— ao^l) 



(D4) 



That d\(t) is regularized "differently" in Eqs. (Dl) and (D2) is a peculiarity of the proper time regularization. In the 
Pauli-Villars regularization all quantities F = {Mn, I, M2{t), J{t), di (t)} are regularized in the same way as 



i^.eg = F{M) - F{Mpv) , (D5) 



PV 
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where it is understood that the corresponding model expressions are first evaluated with the Hamiltonian (26) and 
with the Hamiltonian (26) where M replaced by Mpy, and then finally subtracted according to the prescription (D5). 

Notice that for the constraints (50) and (55) to be satisfied by the numerical results it is of crucial importance 
that M2{t) and Mm as well as J{t) and / are regularized in the same way. This is the case for both regularizations. 
In order to make this apparent also for the proper-time regularization we recall that in this method the regularized 
model expressions for continuum contributions to the nucleon mass (27) and the moment of inertia (30) are given by 

Mjv.cont = ^ {En - £^no)|rcg = ^^c ^ ( i?l(-B„, A) - i?i(£;„,„ A) J 

/eon. . f E = f E (-k»i?3(i?„,i?™, A) . (D6) 

E„<0 m^ii 

Hereby A or Mpv are fixed to reproduce the physical value of the pion decay constant = 93 MeV [29] . For 
M = 350 MeV one has A 649 MeV in the chiral limit, and A = 643 MeV for = 140 MeV. The proper time 
method can, in principle, be applied to any [55]. The Pauli-Villars method can be applied unambiguously in the 
chiral limit, with Mpy = 556 MeV reproducing the experimental value of f-^. However, the method meets difficulties 
in the case of a non-zero current quark mass m [73]. Notice that A and AIpv are of ©(p^^) 600 MeV. 

APPENDIX E: MODEL EXPRESSIONS FOR di(0) AND di'(O) FROM THE EMT 

One consequence of the conservation of the EMT, see Sec. II and Appendices B and C, is that there are two different 
expressions for di(<). In one di(t) is related to p{r) and in the other to s{r). Both are, of course, equivalent. However, 
this is not obvious from the explicit model expressions. 

Here we derive model expressions for di(0) and d[{0) from the EMT in terms of p{r) and s{r) which will be useful 
below for the explicit demonstration that the expressions for di{t) obtained from the EMT and GPDs are equivalent. 

The model expression for the spatial components of the static energy-momentum tensor (8) reads 

= § E C(r)(-7VV+7V^' + (^--0)'/'n(r). (El) 

n, occ ^ ^ 

For the model expressions for the form factor di (t) and its derivative at t = we obtain 

n. OCC ^ ^ n, occ ^ ^ 

dm ^ ^ E (h{/7P, rA\n) = Y: ("|{7°7V, (r'r^ ~ j S^') ]\n) . (E3) 

Ti,occ ^ n, occ ^ ^ ^ ^ 

The first relations in (E2) and (E3) are more practical for a numerical evaluation. They follow for example from 
taking the limit A' ^ in (41) and making use of (16) and hedgehog symmetry. Notice, that (E2) can alternatively 
be derived from (17) and (44) — which is a cross check for the intermediate model expressions. 

In order to find the equivalent second relations in (E2) and (E3) we inspect Eq. (9) for small momentum transfers 
using for T'^ (r) the model expression in (El). This yields 

dm+tldm+oif)^-^^j2(^\{^'^'p'^ (f'f'-^^sA (i + ^+o(t^))|in), (e4) 

n,occ ^ ^ ^ ^ ^ ^ 

from which we read off the desired results. The second relation in (E3) follows also immediately from (13) and (El). 
Comparing (El) with the general relations for p{r) and s{r) with T^^{r), 

p{r) = i T^^r) 5^= , .(r) = ^ (r) {^-^ - U^^^ (E5) 

we recognize that (E2) is just Eq. (17) in the model, while (E3) is the model version of Eq. (68). 
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APPENDIX F: POLYNOMIALITY OF (if" + H''){x,i,t) AT t ^ 

The expression for {H'^ + H'^){x, t), which in the SU(2) version of the CQSM aheady exhausts the sum over quark 
flavours in Eq. (3), was derived, evaluated and discussed in [38]. With the light-like vector n in Eq. (2) chosen to be 
n'^ oc (1, 0, 0, —1) the model expression is given by 

(i/" +i/'^)(.T,e,i) = MnNcJ ^ J2 e*"°("*^"-^")(n| {1 + -f°-f^) exp{-i^f) exp{iAi) exp{-i^f) \n) . (Fl) 



n.occ 



The study of Ref. [38] was supplemented in [40] by the explicit proof that (i?" + H'^){x,£^,t) in the model satisfies 
the polynomiality condition (5) at t = 0. In this Appendix we generalize the proof of Ref. [40] to f 7^ 0. 

The reason why the proof of [40] was restricted to the case t = is connected to the fact that the information on 
5 and t on the RHS of Eq. (Fl) is implicitly encoded in the 3-vector A. In the large-A'^c kinematics A'^ = —2£^Mn 
and t = — A^. By keeping ^ 7^ and continuing analytically the moments of {H^ + H'^){x, t) to the point t = one 
obtains model expressions depending on ^ only. That simplifies the situation considerably [40]. 

However, the proof of [40] can be generalized to i 7^ by using the following remarkable identity 

exp{iAK)=f2i'{2l + l)ji{rV^t)Pi{cos0)Pi(-^^) , (F2) 
1=0 \ V t / 

where r = |x| and cos9 ~ a;'^/|x| and ji{z) and Pi{z) denote respectively Bessel functions and Legendre polynomials. 
As a byproduct we remark that by inserting the identity (F2) in (Fl) we see explicitly that (ff" + -ff'')(.T, t) IS a 
fimction only of f and t (and x, of course) but not of the 3-vector A. 

The identity (F2) can be derived as follows. We rewrite the 3-vector A as A = (sin (5 cos a, sin (5 sin a, cos 5) -\/^ 
where cos (5 = -~2£^M]y/^/^ and a £ [0; 27r] describes the orientation of A in the plane perpendicular to the spatial 
direction of the light-cone. Then, by using spherical coordinates for x, we obtain 

00 

exp(jAx) = ^i'(2/ + l)j((r\/^)^';(cos6'cos(5 + sin6'sin(5cos(a - 0)) . (F3) 

(=0 

Next we rewrite the Legendre polynomials using the addition theorem as 

Pi (cos 61 cos 5 + sin 9 sin 6 cos(a ~ cj))) ^ Pi (cos 9)Pi (cos (5) + 2 ^ 77"^ ^)^/ ("^^^ cos(/(a - 0)) , (F4) 

The physical process described by the GPD does not depend on the orientation of A in the transverse plane. This is 
reflected by the fact that the matrix elements in (Fl) are independent of a. (Changes of a can be compensated by 
appropriate hedgehog rotations.) Thus, we can eliminate the artificial dependence on a — for example by averaging 
the matrix elements in (Fl) over a. The latter is equivalent to taking the average over a in (F4). Inserting the result 
into (F3) gives finally the identity (F2). 

Notice that for e ^ the product {21 + 1) ji{Ae) Pi{B/e) -> {ABY/U + ©(e^). Thus, taking the limit t while 
keeping ^ 7^ in (F2) yields 

Jim^^ exp(*Ax)-f;^^^^^^j^^'p,(cos0) , (F5) 

reproducing Eq. (20) of Ref. [40] which was derived independently in a different way. 

The proof of polynomiality at i = given in Ref. [40] is generalized to any t 7^ by replacing Eq. (20) in [40] by 
our more general identity (F2) and repeating literally the steps done in Eqs. (21-28) in [40]. This yields 

/d.."-(i/" + i7^)(.,^,t) = ^ X: ^\2i + m('-^ 

fe 

X 



E ■ ("i(7V)' {pyji{mv^t)Piicos9) ip')''-^n) , (f6) 
3=0 ^ ^ 

and explicitly demonstrates that the model expression (Fl) does satisfy the polynomiality property for any t 7^ 0. 
Note that this includes also positive t allowing to study in principle the region of time- like momentum transfers, where 
GPDs "become" nuclcon-antinucleon distribution amplitudes [83, 84]. 
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APPENDIX G: M2{t) AND di{t) FROM GPDS 

111 this and the foUowing Appendix we introduce labels to distinguish the model expressions for form factors derived 
from GPDs and EMT — with the aim to prove finally that they are equivalent. From Eq. (F6) we read off the model 
expressions for the second Mellin moment of (iJ" + H'^){x,£,,t), see Eq. (3), 



M2WIGPD = ^EHi?r,Jo(|x|V^)|n) + ^^ 



n.occ 



{7VJ5^Jo(|A|^/^) + ^ j2(|x|y=i)P2(cos^)} 



rflWlGPD = 75iV,MAr ^ {j0^^f,j,(\S,\^t)P,{coBe)} \n) . (Gl) 

n,occ 

By comparing the expressions for -/\/2(^)|gpd and (ii(t)|GPD in (Gl), and by exploring hedgehog symmetry we observe 

^ n,occ ^ 

{n\EM^\V^t)\n) = ^ E ("l^-e'^*!'^) ^ iM^it) - * di(i)]EMT 

Mm ^ — ^ Mm ^ — ^ '^Mj.r 

n,occ n.occ -'^ 

^ E {n\{l'l'f:Joi\i\V^t)}\n) = ^ E ("l{7"7P,e^^*}l") ^ y^^iWemt - (G2) 

n.nrr n.nr.r. ^ 



Thus, Eq. (Gl) states that 



[M2{t) - TTTTT rfl WIgPD - [M,it) - — dl(t)]EMT . (G3) 



This means that the model expressions from GPDs and from EMT are equivalent for this particular linear combination 
of M2{t) and di{t). In order to prove that this is true also for the separate form factors we have show explicitly that 
e.g. one can derive the same model expression for (ii(i:)GPD from the EMT. One way to do this is to demonstrate that 
'^i(Ogpd satisfies the differential equation (10) derived from the EMT with appropriate boundary conditions. 

For that let us first remove the preference of the 3-axis from the expression for (ii(t)|GPD in (Gl) which is due to 
arbitrarily choosing the light-like vector (x (1,0,0,-1), see above. We remove this arbitrariness by averaging the 
expression (Gl) over directions. This yields 



and using the identity 



4^ d ^ J2(|x|^/=I) _ jo(|x|V^) 

3 dt 15 dt^ J |X|2(-t) 15 ^ ' 



we see that fii(t)GPD satisfies the differential equation (10) 

[d,{t) + jd[{t) + ^d['{t)]aPU = EH {7°7V, U^i'^l\S.\W\jo{\S.\V^t)\\n) 

= yd^re— ^r^(r) j^rV^ - ^ (5''^^ . (G6) 

Next, using (G4) and expanding the expression for (ii(i)GPD in (Gl) we obtain 

rfi(o)iGPD = -^^Eh{7°7¥, Ui'^lm's'^ 

n.occ ^ ^ 



rf{(0)|GPD = - ^ E (-1 {7°7V, Ui' - I m'S'') |xp| In) , (G7) 

n,occ ^ ^ ^ 

which coincides with the expressions obtained from the EMT in Eqs. (E2) and (E3). This, completes the proof that 
one obtains the same model expressions for the form factors di{t) and M2{t) from GPDs and from the EMT. 
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APPENDIX H: J{t) FROM GPDS 

In this Appendix we show that the model expression for J(t)EMT obtained from the EMT in Eq. (42) coincides 
with the model expression for J(f)GPD which results from GPDs [3] by adding up the sum rules (3) and (4). 

The model expression for (77" + Ef^ + + E'^){x, ^, t), which wc shall refer to as Em{x, t) for brevity, reads [42] 



E 



-iz"E„ 



E,,, — Ea 



X {rn\T^\i) {i\ (1 + 7°7^) cxp(-izV/2) cxp(iAX) cxp(-izV/2) \m) 

-0 E e-'"°^'"(m|T''(l + 7V) exp(-izV/2) exp(zAX) exp(-izV/2) |m) 



-iz 2^ e 

m.occ 



(HI) 



Integrating the xEi\j{x, S,,t) over x yields, after substituting x y ^ xMjv and extending the limits of y-integration 
to the entire y-axis which is justified in the large A^c-liniit due to AIn = 0{Nc), the following result 



da; xEm{x, t) = — 



E ^™ - E ^1 

m,occ m.an 

J. all j.occ 



1 



E — E- 



Hr'\j) (j|(l + 7V) exp(zAX)lm) 



Em — E,i 



MA]) (1 +7V){P^cxp(^AX)} \m) 



{m\T^ (1 + 7%^) cxp(iAX) \m) 



(H2) 



Consider the unitary transformation G5 = 7'^7^r^ in the notation of [85] with the properties G^^^G^^ ~ {l^)'^ 
GsT^Gg"^ = — (r°)-^. It transforms in coordinate space the Hamiltonian (26) and its eigenstates as G^HG^^ — 
and G,5$,i(x) = $* (x). Making use of this transformation wc obtain the identities 



one obtains 



(j|rV) = 

01(1+ 7V) exp(*AX)|m) = 

{j\T^ (1 + 7%^) exp(iAX) |m) = 

01(l+7V){p',exp(*AX)}|TO) = 



(m|(-r^)|j) 

(m|(l-7V) exp(zAX)lj) 
{m\T^ (-1 + 7°73) exp(iAX) \j) 
(m|(l-7V) {-p3,exp(*AX)}|j) 



1 

j dx xEM{x,i,t) ^ E 



(H3) 



(H4) 



with 



Em — Ej 
^ ^abc 

3 Em — Ei 



Kfn\T \]) 01 



'b) 01 



f-l\Em + S,) exp(i AX) + {p\ exp(^AX)} 



l^'-I^Em + Sj) exp(i AX) + {p^ exp(zAX)} 



(H5) 



where in the second step we made use of the hedgehog symmetry. The final step necessary to recognize that (114, H5) 
coincide with the expression (42) is the following. Notice that in the expression in Eqs. (H4, H5) there is no more 
reference to the 3-direction which was picked out by arbitrarily choosing the space-direction of light-cone vector 
along the 3-axis. Therefore it is justified to identify A^ = |A^ = — 1 1, which finally leads us to the result 



2J(t)GPD = j d.T xEM{x,i,t) = 2J(i)EMT 



(H6) 
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APPENDIX I: CHIRAL PROPERTIES OF THE FORM FACTORS 



111 this Appendix we study the chiral properties of the EMT form factors, and start with di{t) because here we face 
simpler expressions and the method is more easily explained. We start from the expression (Gl) for di{t), and rewrite 
its continuum contribution in terms of the Feynman propagator, see App. A of Ref. [32], and expand it in gradients 
of the f/-field. In leading order of this expansion we obtain 

rfiWcont = ^^4r^ yd^x j2(rx/^)P2(cos0)trF[V3t/][V3t/t] + ... (II) 

where the dots denote terms containing three or more gradients of the ?7-ficld. Taking i ^ in (II), and using 
i'(^) = (2iii)V. + C>{z^+^) for z < 1, we recover the result for d^""* in Eq. (44) of Ref. [40]."* 

The result in Eq. (II) can be used to study the chiral properties oi di{t). For that the leading large-r (long distance) 
behaviour of the integrand in (II) plays the crucial role. Therefore for this purpose it is legitimate to neglect both, the 
discrete level contribution which has exponential fall-off at large r, and higher order terms in the gradient expansion 
which have additional power suppression with respect to the leading order term in (II). 

In Eq. (Gl) we have chosen as starting point di{t) is given in terms of s{r), see App. G. Therefore — after taking 
the flavour-trace trp in (II), integrating out the angular dependence, restoring the integral over the full solid angle, 
and comparing to Eq. (17) — we read off the expression for s(r), 

di(Oco„t = ^ jd?^n{rV^t)s{T) with s{r) = fl{p'{rf-^^^^^^+... (12) 

Making use of the long-distance behaviour of the profile (29) we obtain for the large-r behaviour of s{r) in the chiral 
limit the result quoted in Eq. (62). The large-r behaviour of p(r) quoted there follows from (15). 

In order to derive the leading non-analytic contributions to di{t) we choose to work with the following analytic 
form of the soliton profile 



P{r) = —2 arctan 



2 



H 

— (l-f m^r) exp(-m^r) 



(13) 



This profile was demonstrated to be a good approximation to the true self-consistent profile [26]. Since it does not 
correspond to the true minimum of the soliton energy, e.g., the approximate result for the pressure obtained with this 
profile does not satisfy the stability condition (16). However, all that matters for our purposes is that it exhibits the 
correct chiral behaviour, see (29). The "soliton radius" R in (13) is connected to the constant A in (29) by ^ = 2R^. 
Let us focus on the chiral expansion of rfi(i)cont at zero momentum transfer. We obtain 

oo 

^1 = -Y-f'^^NR'Gim^R) with G{a) = jdz (^4 + (l + az)^^-^^^)^ ^^^^ 



The zeroth order in the Taylor series of the function G[a) around a = is G(0) = J^^z 12 z^ / (z^ -I- 1)^ = 97r/-\/8. 

To find the linear term in the Taylor expansion of G{a) we proceed as follows, c.f. App. B of [31] for a similar 
calculation. We consider G'{a) at a 7^ 0, substitute z ^ y = az, and consider then the limit a ^ 0. Notice that, had 
we decided to evaluate (14) in a finite volume, let us say in a spherical box of radius D, then the upper limit of the 
y- integration would be DttIt^. Thus, before taking the chiral limit it is crucial to take first the infinite volume limit 
-D — > cx) [31], as these two limits do not commute. This yields G'{a)\a=o ~ ~8 J^dy (—1 + y + 2y^ + y^) e'^^ = —5. 
Thus, we obtain 

G{a) = G{0) — 5a + higher order terms in a. (15) 



* In [40] in the leading order in gradient expansion it was estimated i^j'^g^!^^"* ~ —9.46 to be compared to rfj°"xact ~ —8.34 which is the 
exact result obtained here in proper time regularization, i.e. also in this case the gradient expansion provides a useful estimate for the 
continuum contribution of a quantity. Let us take here the opportunity to correct an error in [40] . The level contribution to di is not 
zero — contrary to the claim in Eq. (43) of [40], Instead, it is = 4.88 with the self-consistent proper-time profile. In fact, the level 
contribution to di (and to the pressure) cannot be zero. It plays a crucial role in establishing the stability of the soliton, see Sec. VII. 
We stress that reliable results which satisfy the stability condition and all other requirements can only be obtained from evaluating 
numerically the exact model expressions, with the correct self-consistent profile, see Sec. VII. 
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Inserting this result in (14) we reproduce the chiral Hmit result in Eq. (46) of [40] which, however, provides just one 
contribution to the chiral limit value of di. Other contributions are of importance, see Footnote 4. 

The situation is different for the linear m7r-correction to c?i. As explained above, here the expansion (15) provides, in 
fact, the correct leading non-analytic term in the chiral expansion of the full expression for di in the model. Combining 
(14) and (15) and eliminating R = y/ A/2 in favour of gA and /tt according to (29) we obtain Eq. (65). The result for 
d[{0) in (70) is obtained similarly. 

To derive the small t expansion of di{t) in the chiral limit in Eq. (69) one can integrate (12) exactly. The result is 
a bulky and not illuminating expression which we do not quote here. Expanding it for small t yields (69). 

Let us turn to the discussion of the form factor M2{t). At t = we have for all m^r exactly M2(0) = 1, see Eq. (50). 
The slope of M2{t) at zero momentum transfer, however, has a non-trivial chiral expansion. Due to (66) we need for 
that the chiral expansion of di (t) , see above, and that of the mean square radius of the energy density. In the above 
described way we obtain from the expression (52) for the energy density in the leading order gradient expansion the 
following contribution 



(4)partl = (?|)-f (4./.^i?4)^-f- 



(16) 



which is, however, not yet the complete result for the following reason. The chiral expansion of the exact model 
expression for the energy density (43) in the gradient expansion contains — in addition to Eq. (52) — also the term 



tTFi2 - U - U^) ^ mlf^ll - cosPir) 



(17) 



This term is of "zeroth order" in the gradient expansion. It arises from the current quark mass term in (25) and is 
related to the nucleon-pion sigma-term in the gradient expansion [31]. The explicit appearance of the current quark 
mass has been eliminated in (17) in favour of by means of the Gell-Mann~Oakes~Renner relation which is valid in 
the model. The mj-term (17) vanishes in the chiral limit, and was therefore not displayed in (52), but it contributes 
in linear TO7r-order to (16). This may not be obvious at a first glance, however, from (17) we obtain 



' E/par 



oo 



z«{z^ - 2z - l)e-2^ - 4a4z4(l + 



{z* + a^{l + zfe-'^')'^ 



a=0 



(18) 



and adding up (16) and (18) we obtain 



(4) 



{rl) 



9(4./^i^^)^ + 



(19) 



It is not necessary to consider chiral corrections due to the nucleon mass Mn = Mn +Bm^ + . . . because they 
contribute only at higher orders. (Notice that the CQSM also consistently describes the chiral expansion of [31].) 
Ehminating R = y/A/2 in (19) by means of (29) yields finally the results for (r|,) in (53) and for M2(0) in (66). 

Finally we comment on the form factor J{t). Also in this case the normalization is trivial, since J(0) = i, but 
e.g. the chiral expansion of J'(0) = ^ (rj) is of interest. Here we restrict ourselves to the mere observation that (rj) 
diverges in the chiral limit. 

In contrast to M2{t) and di(t), which are non-zero in the leading order of the large- TVc expansion, J{t) arises from 
1/A^c ("rotational") corrections. For such quantities the non-commutativity of the large- A''c and chiral limit may have 
more drastic implications [48]. For example, in the slowly rotating soliton approach (as realized e.g. in the Skyrme 
model in Ref. [86]) the isovector electric mean square radius diverges as l/mj^ in the chiral limit — in contrast to 
Inm^r at finite Nc- For (rj) the situation is completely analog — as a study in the Skyrme model reveals [87]. 
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